Real-world Quantum Sensors: Evaluating Resources for Precision Measurement 
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Quantum physics holds the promise of enabling certain tasks with better performance than pos- 
sible when only classical resources are employed. The quantum phenomena present in many experi- 
ments signify nonclassical behavior, but do not always imply superior performance. Quantifying the 
enhancement achieved from quantum behavior requires careful analysis of the resources involved. 
We analyze the specific case of parameter estimation using an optical interferometer, where increased 
precision can be achieved using quantum probe states. Common performance measures are exam- 
ined and it is shown that some overestimate the improvement. For the simplest experimental case we 
compare the different measures and show this overestimate explicitly. We give the preferred analysis 
of real-world experiments and calculate benchmark values for experimental parameters necessary to 
realize a precision enhancement. 



Quantum-enhanced technologies such as quantum 
computing, cryptography and metrology utilize nonclas- 
sical behavior of quantum systems to surpass the per- 
formance of their classical analogs. Often signatures of 
quantum behavior such as violation of a Bell inequality, 
photon antibunching, or sub-Poissonian number statis- 
tics are observed in an experiment. In some cases, their 
presence may be sufficient to ensure superior performance 
of the quantum protocol over the classical. There are sit- 
uations, however, in which quantum behavior does not 
necessarily imply performance beyond classical limits and 
further analysis is necessary to ensure improvement. 

An important class of quantum-enhanced technologies 
involves sensing, in which quantum probes are used to 
estimate a parameter of a physical system with preci- 
sion beyond what is possible with classical resources [l[. 
Analysis of the limits on precision, something which is of 
both practical and fundamental interest, requires careful 
accounting of the necessary resources to compare classi- 
cal and quantum protocols. Entangled states of multiple 
particles are a critical component in enabling quantum 
enhancement in many situations ^ , as are detectors that 
extract maximum information about the parameter of in- 
terest [3|, what we call optimal measurements. Calculat- 
ing the precision that can be achieved for fixed resources 
provides a basis for comparing quantum and classical pro- 
tocols. Recent work has shown that practical imperfec- 
tions such as loss, decoherence, state-preparation and de- 
tector inefficiency decrease the enhancement of quantum 
metrology protocols [^-|6|. Therefore poor enumeration, 
which distorts the resources used and hence the enhance- 
ment obtained, is problematic. 

Several proof-of-principle experiments aimed at 
quantum-enhanced parameter estimation have been per- 
formed to date [Tj]. Many utilize nonclassical signatures 
in detection outcomes to indicate improved precision. 
However, as we will show, this is insufficient to demon- 
strate increased precision beyond the classical limit and 
further information must be included to accurately com- 
pare quantum and classical strategies. 



In this Letter, we examine a quantum sensor based 
on optical interferometry. We evaluate the resources, in 
terms of number of trials and photons, to reach a given 
precision when estimating the phase between two opti- 
cal modes. Limits to this precision are analysed using 
Fisher information and quantum Fisher information and 
compared with other common measures of enhancement 
such as super-sensitivity [8|. When experimental imper- 
fections are present, we show that quantum enhancement 
can be overestimated by improper resource accounting. 
In particular, a common error is to neglect the channel 
transmission, source, and detection efficiencies. To ex- 
perimentally demonstrate this, we construct the simplest 
quantum states which ideally lead to improved precision 
and examine the behavior of different enhancement mea- 
sures. We show that improvement in state preparation 
and detection efficiencies are necessary to unambiguously 
demonstrate quantum enhancement in phase estimation, 
and set quantitative benchmarks for these. 

To estimate the phase (/) in a two-mode interferometer, 
a state p is launched into it, evolves into the state p((/)) 
due to interaction with the phase-shifting element, and 
a measurement is performed at the sensor output. The 
phase (j) is estimated from the outcomes of v repeated ex- 
perimental trials. We define the average number of pho- 
tons at the interferometer input, N = (n) , and the total 
number of experimental trials, v, necessary to achieve 
a given level of phase precision, Acj), as the resources. 
For fixed resources, we then compare the precision of the 
quantum and classical approaches and the approach with 
the smallest phase uncertainty is thus the better strategy. 

Two effects attributed to quantum behavior in inter- 
ferometry are phase super-resolution and phase super- 
sensitivity. Phase super-resolution, the sinusoidal varia- 
tion of an A^-fold detection signal as a function of inter- 
ferometer phase cf> with an iV-times increase in rate of os- 
cillation, can be observed using only classical input states 
of light and projective measurements [8|. It has recently 
been shown that the visibility of such super-resolving 
fringes can be used to distinguish between quantum and 



classical input states of a Mach-Zchnder interferometer 
[9|. This signature of quantum behavior does not, how- 
ever, quantify improved performance beyond classical in- 
terferometry. 

Phase super-sensitivity, a commonly employed mea- 
sure of performance, is defined as reduced phase un- 
certainty compared to that possible with classical re- 
sources [8| . The model employed in Ref. [8| incorporates 
experimental imperfection through 'efficiency' and visi- 
bility parameters ryp and V , by means of a phenomeno- 
logical model of an A^d-fold detected coincidence signal. 
Here, rjp is the proportion of the input state p that can 
lead to an iVd-fold detection event. The visibility is re- 
quired to satisfy i^pV^Nd > 1 in order for the measure- 
ment to be regarded as 'super-sensitive'. However, this 
assumes unit sensor transmission and perfectly efficient 
detectors. It includes only imperfect state preparation; 
with N^ = N, the number of input photons. In fact, the 
most common experimental situation is that A^d < N 
so that the resources consumed are significantly under- 
estimated. Claims of precision beyond the classical limit 
that utilize this measure must therefore be interpreted 
carefully. 

In general, the precision of the phase estimate is 
limited by the Cramer- Rao bound (CRB) [3|, Ai/) > 
l/yiZF(0), where F(</.) = E,{'^/pM)\dp,{^)/dc^\\ is 
the Fisher information (FI). The probability Pj{(j)) cor- 
responds to an outcome j of a measurement. F{(l>) is 
bounded from above by the quantum Fisher information 
(QFI), Fq, found by maximizing F{(j)) over all physical 
measurements. Hence the quantum Cramer-Rao bound 
(QCRB), A(/) > 1/^:/Fq((/)), depends only on the in- 
put state and channel, independent of detector configu- 
ration [3[. 

In the absence of imperfections, the QCRB of N un- 
correlated photons leads to the standard quantum limit 
A(j) > A(/)sQL = l/yJvN, equivalent to the precision at- 
tained with a coherent state of unknown phase with av- 
erage photon number |ap = A^. When quantum input 
states are employed, the best attainable precision is given 
by the so-called Heisenberg limit /S.(f) > A0hl = 1/v^-^ 
which is achieved by the widely studied 'iVOOiV state' [7|- 
|lO|. Note that to reach this precision, an appropriate 
optimal detection scheme with unit efficiency must be 
used. However, if the sensor has finite transmissivity 77, 
assumed equal in both interferometer arms, the classical 
strategy is bounded by the standard interferometric limit 
(SIL) A(t> > A(/)siL = l/y/T^ H. A^OOA^ states are ex- 
tremely susceptible to loss; removal of even one photon 
results in a phase-insensitive state. The QFI for A^OOA^ 
states through the same channel is thus scaled by the 
probability that all photons are transmitted, r; , giving 



This measurement set has three possible outcomes: j ~ 
± (detection of |^±)) and j = (otherwise). In practice, 
it is difficult to realize single-mode sensors and this may 
lead to non-ideal interference due to unmeasured distin- 
guishing information. Measurement outcomes therefore 
typically have probabilities p±[(j)) = /[I ± V cos{N(j))]/2 
and po[(l)) = 1 — p+(0) — V-jct)), leading to a CRB 
Acj, > A</.„,i„ = l/V^NV [ni. Here, / = Vpivm)"" , 
where Tyd is detector efficiency, and V is the interference 
visibility accounting for multimode states and detectors. 
The state is said to exhibit super-resolution if p±((/>) os- 
cillates N times the applied phase. To surpass the SIL 
requires Acjj^i^ < A0giL, leading to a threshold visibility 



T4h = Vv/fN, 



(1) 



which is equivalent to that in Ref. [8| with 77 = yyd = 1 
and N = A^d- A state is said to exhibit super-sensitivity 
ii V > Vth, and the classical precision limit is surpassed 
when this occurs. 

The FI and QFI analyses give a straightforward 
method to quantify the resources consumed in an ex- 
periment and compare quantum with classical strategies. 
The QFI reveals whether or not the classical limit could 
be beaten in principle using a particular state and chan- 
nel, while the FI reveals whether or not it can be beaten 
with the addition of a particular measurement. To ascer- 
tain the best performance of a given quantum strategy, 
one must therefore know both the density matrix p of the 
input state, giving TV, FI, and QFI, and when it is pre- 
pared, giving v, as well as the device transmission and de- 
tector efficiency. In our case, p will be obtained through 
state tomography and i^ by heralded state preparation. 

A common practice is to 'post-select' on particular 
measurement outcomes and neglect the occurence of oth- 
ers, including when nothing is detected at the output. 
This is equivalent to setting ?/p, 77, r]^ to 1 and A^ = A^d- 
This neglects both v and the true A^, significantly un- 
derestimating both the number of trials (which grows 
exponentially in N [ll|) and information from additional 
measurement outcomes. 

To experimentally examine realistic quantum sensors, 
we implement the simplest case, employing two-particle 
states, in which increased precision can be shown. We use 
an optical sensor design based on a heralded two-photon 
Holland-Burnett (HB) state [il|,[il. HB(A^ ::= 2k) states 



a QCRB of A(j) > 1/y/mfN. 

For NOON states, an optimal measurement is projec- 
tion onto |^±) = (|A^,0) ± \0,N))/V2, where \m,n) de- 
notes m (n) photons in each mode of the interferometer. 



are prepared by launching k photons in both input ports 
of a 50:50 beam splitter. The HB(2) state provides the 
least stringent constraints on sensor transmissivity and 
detection efficiency [J,l5|. Furthermore, the optimal mea- 
surement for phase estimation using HB states can be 
implemented with realistic number-resolving detection 
[5|. The QCRB of a HB state in the lossless case is 
A(j) > l/^uN{N/2 + 1) which, for large A^, has Heisen- 
berg scaling that differs only by a constant factor. With 
loss this state retains phase sensitivity and precision en- 
hancement far better than A^OOA^ states [13|. Indeed, 



the QCRB for the HB state approximately foUows the 
optimuni value for a wide range of ry, rj^ and iV [J, |5|. 

The scheme used to generate heralded HB states is 
shown in Fig. [TJa). A polarization-based iV-photon HB 
state is obtained by combining two orthogonally polar- 
ized k = 7V/2-photon Fock states at a polarizing beam 
splitter 111|. The required Fock states are generated by 
two parametric downconverters (PDC), where k photons 
in one mode of each PDC are heralded by detection of 
k photons in the other mode. The PDCs are designed 
to produce spectrally disentangled modes [14| so that 
heralding does not compromise the state purity, which 
is crucial for high-visibility interferometry. Since spec- 
tral filtering is not required to increase the state pu- 
rity, the preparation efficiency is determined solely by 
the system loss. In principle this source, with the use of 
photon- number- resolving heralding detectors |15| , gener- 
ates ideal HB states of arbitrary N. 
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purity and distinguishability, and an upper bound for the 
multi-photon fringe visibility in Eq. ([1]) . The residual im- 
purity is partly intrinsic [IJ] and partly due to imperfect 
compensation of the optical fiber birefringence. 

We reconstruct the full heralded state by performing 
state tomography u sing the experimental arrangement 
shown in Fig. [TJc) [ll|, ll^. Population in lower pho- 
ton number subspaces arises from loss so that coherences 
between subspaces of different photon number can be ne- 
glected. Figure [Slja) shows the reconstructed density ma- 
trix. The heralded state has populations of 0.686, 0.277, 
and 0.037 in the zero-, one-, and two-photon subspaces 
giving an average photon number of 0.35. This state has 
an overlap of 0.031 with the ideal 2-photon HB state, 
while the re-normalized two photon subspace has an over- 
lap of 0.85. 

The heralded HB states are launched into a polariza- 
tion Mach-Zchnder interferometer, PMZI in Fig. [TJb), 
which introduces a controllable relative phase between 
the two modes. Figure [U^b) shows the detection count- 
rate phase dependence for heralded single-photon states 
(dashed) and heralded HB(2) states (solid). Theoreti- 
cal fits show visibilities of 90.9 ± 0.5% and 80.6 ± 1.5% 
(76.4 ± 1.4% raw) respectively. The latter is consistent 
with the measured state purity. 



FIG. 1. (color online) Setup for generation and characteriza- 
tion of heralded HB states, (a) State generation based on two 
parametric downconverters (PDC). Half- (H) and a quarter- 
wave plates (Q) after each polarizing beam splitter (PBS) 
adjust the polarizations to combine photons at a fiber polar- 
izing beam splitter (FPBS). Coincidence detection between 
heralding avalanche photodiodes (Herald APDs) signals HB 
state preparation at the output of the FPBS. (b) PMZI ap- 
plies phase (j) between ±45° polarizations. Coincidence de- 
tection between the APDs implement optimal measurement 
for HB(2) state, (c) Tomographic measurements are imple- 
mented using a quarter- (Q) and half-wave plate (H) followed 
by a PBS. Outputs are coupled into single-mode fiber (FC) 
with the reflected mode split by a 50:50 fiber beam splitter 
(FBS) to allow partial number resolution. 

The quality of the heralded Fock states used to gen- 
erate the HB(2) state was tested by means of Hong- 
Ou-Mandel interference between heralded single photons 
(Fig. Hi). The theoretical fit gives a visibility of 90 ± 3% 
(80 ±3% raw), setting a lower bound on the input photon 



FIG. 2. (color online) Quantum intereference of heralded pho- 
tonic states, (a) Hong-Ou-Mandel interference between two 
heralded single photons as a function of delay, r. (b) Set- 
ting T — and scanning the PMZI phase shows two-photon 
(solid) and single-photon (dashed) interference. Error bars 
are derived from Poissonian statistics and contributions from 
multiple PDC pair emissions are removed. 

For N = 2, HB and A^OOA^ states are equivalent and 
Eq. (ID) can be applied. If we use a post-selection analy- 
sis, only data resulting in 2-fold outcomes are recorded, 
■q and / are set to 1 and N = 2. This yields a threshold 
visibility Vth = l/\/2. The measured fringe visibility sig- 
nificantly exceeds this bound and so by this analysis we 
achieve super-sensitivity, beating the SQL. If we renor- 
malise the post-selected photon number subspace and 
use r/p = 0.85 while still neglecting 77 and yyd, we obtain 
Vth = 0.77. Our measured fringe visibility also exceeds 
this threshold, again demonstrating super-sensitivity and 
beating the SQL. Since the whole state including lower 
photon number subspaces is known however, the correct 
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FIG. 3. (color online) Reconstructed heralded quantum state, 
(a) The absolute value of the density matrix including vac- 
uum, one-, and two-photon subspaces. (b) An enlargement 
of the two-photon subspace. In all density matrices, \m, n) 
denotes m diagonally (-f45°) polarized photons and n anti- 
diagonally (—45°) polarized photons. Bars are colored ac- 
cording to their magnitude. 

77p, given by the overlap of the reconstructed state with 
the HB(2) state, can be calculated as Typ = 0.031. We 
can therefore apply the bound in Ref. [8[, leading to 



Vth = l/v'7p-^d = 4.02. This is an unphysical visibility 
and so, by this analysis, we cannot beat the SQL. The 
above threshold visibilities are based on post-selection 
and renornialisation or neglect both channel transmis- 
sivity and detector efficiencies. They therefore signifi- 
cantly underestimate the resources required to reach a 
given precision. Heralding a characterised state affords 
a complete reckoning of the resources used by counting 
and all states input into the interferometer. / and rj can 
then be estimated experimentally, allowing Vth to be cal- 
culated according to Eq. ([1]). / is determined from the 
ratio of four-fold coincidences to heralding events, giv- 
ing / = 0.0047 ± 0.0001 and yyyyd = y/TTrh> = 0-39. For 
our detector efficiencies of r]p = 0.45%, rj = 0.87 which 
results in Vth = x/v/^f ~ 9.57, an even higher thresh- 
old on visibility, showing that we cannot beat the SIL. 
This explicitly shows that to unequivocally demonstrate 
performance beyond the classical limit, all system ineffi- 
ciencies must be characterised. 

Equation ([1]) and the final analysis above arc based on 
the FI and arc valid only for iVOOiV states. In general, 
however, to assess the value of a state for metrology the 
QFI can be calculated and compared to the SIL. For our 
state, Fq = 0.079, giving a QCRB A0 > 3.56/v^. For a 
classical state with the same average photon number as 
the reconstructed state, Acfisih = l.Sl/y^. This shows 
again that, despite the implications of the post selected 
and renormalizcd analyses above, the classical limits to 
precision can never be surpassed with this state and ar- 
bitrary measurements. 

Going beyond the classical performance limit with 
realistic quantum sensors requires stringent bounds on 
throughput, probe-state preparation and detection effi- 



ciencies. For heralded NOON states, / = ViAvid) must 
satisfy / > rj/NV"^. For the ideal two-photon situa- 
tion iy = 1) this implies rjprjrj^ > 1/2. This bench- 
mark for heralded two-photon states provides a challeng- 
ing goal that must be achieved to surpass the classical 
limit. Even for perfect detectors (ryd = 1), this would 
require ijpr] > 1/2. Conversely, perfect preparation and 
transmission {r/pr] = 1) implies that any detectors used 
must have rjd > l/\/2. To our knowledge, the above cri- 
teria have not yet been demonstrated together in a single 
experiment. 

Quantum photonics using feasible laboratory technol- 
ogy has the potential to surpass the performance of clas- 
sical techniques. Nonclassical behavior may be present in 
an experiment, but does not necessarily imply improved 
performance beyond classical limits. Here we have laid 
out an explicit approach to compare experimental results 
for precision metrology, which requires complete knowl- 
edge of the input state and number of experimental trials. 
We have demonstrated a scalable method to generate her- 
alded entangled states of light for precision phase estima- 
tion, a key step to demonstrating a real-world improve- 
ment over the SQL. By fully characterizing this source 
we have calculated the ultimate precision that can be 
obtained and in doing so we have put bounds on the 
sensor parameters required for quantum enhanced phase 
estimation to be viable. 
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Fisher information and phase uncertainty 

The amount of information about the phase 4> between 
two modes in an interferometer that can be extracted by 
a given experimental configuration is quantified by the 
Fisher information 



j^w = E 



1 



pjW 



dpj{(f>) 



(2) 



Here the probabilities Pj{4i) correspond to a particular 
outcome j of a measurement described by a set of opera- 
tors {Hj} when the actual parameter to be estimated has 
a value (j). The precision with which (j) can be estimated 
is determined by the Cramer-Rao bound (CRB) jl7| . 



A0 > 1/^lJFi^, 



(3) 



which can be achieved for a large number of trials v and 
an unbaised estimator, such as a maximum likelihood 
method. The CRB is bounded from below by the quan- 
tum Cramer- Rao bound (QCRB), which is independent 
of the measurement and depends only on the state used, 



A0 > l/V^^Fi^ > 1/^vFq{4>), 



(4) 



where Fn (S) = Fq {p{4>)) is the quantum Fisher infor- 



mation 



rrnigt- Note that in general both the FI and 
the QFI depend on the value of the parameter (f>. Only 
in certain special cases are they independent of the pa- 
rameter. 

In a perfect experiment, an iV-fold detection event cor- 
responds to projection onto an iV-photon A^OOA^ state, 
l*±) = (1^,0) ± |0,iV))/V2, where \m,n) denotes m 
(n) photons in mode a (6). This measurement set has 
three possible outcomes: j = ± (detection of |^±)) and 
j = (otherwise). The probabilities associated with 
these measurement outcomes for a NOON state input are 
P±{4') = [1 i co&{N (j))] / 2 . In practice, it is difficult to 
realize single-mode sensors and this can lead to degra- 
dation of interference due to unmeasured distinguishing 
information of the input state. Measurement outcomes 
therefore typically have probabilities 



p±(0) = /[l±ycos(7V0)]/2, 



(5) 



and po(0) = 1 — P+(0) — P-{<f>), where / and V both 
depend upon the presence of distinguishing information 
in the input state. Here / is related to the probability 
that an input state leads to an iV-fold detection event, 
and V is the fringe visibility. The Fisher information for 
this configuration is given by 

sin^ {N(j)) 



F{cj)) - fN'V 



1 



(6) 



y2cos2(iV0)' 

where we have made use of the first derivatives of p± {(j)) 
with respect to (j), 



^ = ^Wsi„W,. 



(7) 



Note that the denominator in Eq. ^ is greater than or 
equal to s\t?{N(J)) for < 1^ < 1, which sets a bound on 
the Fisher information, 



F{(f) < fV^N^ 



(8) 



The phase uncertainty that can be achieved with ly trials 
is then bounded by the Cramer-Rao bound 



(9) 



A0 > l/i^iyfNV). 



Comparison with the precision that can be achieved in 
a classical approach with N particles, i.e. the standard 
quantum limit (SQL), 



A0SQL = 1/\A^, 



(10) 



implies that in order to demonstrate precision beyond the 
SQL, 

(11) 



Effects of loss and detection inefficiencies on 
measurement precision 

We consider the case of a single spectral-mode sen- 
sor, since this enables / to be calculated straightfor- 
wardly and allows comparison with previous work involv- 
ing super-resolution and super-sensitivity criteria. For 
the case of no loss and perfect detectors, we define f = r/p, 
the production efficiency of the single-mode A^OOA^ state. 
When there is balanced loss in the interferometer, char- 
acterized by the transmissivity rj and detector efficiency 
77d , the probability that an TV-photon A^ 00 A^ state is suc- 
cessfully detected is decreased by a factor (?7?7d)^, giv- 
ing / = rjpirjrjd)'^ ■ In the classical case, the proba- 
bility that a coherent state is transmitted is simply i] 
so that the classical phase uncertainty is bounded by 
A0SIL = 1/ y/vrjN y. This bound is the best precision 
that can be achieved by using the same channel with a 
classical input state and perfect detectors. 

/ can be estimated from the ratio of successful 
A^OOA^ detection events to herald events, giving rjrjd = 
'{f f /rjp. To surpass the performance allowed classically, 
i.e. A(/)siL > 1/(aA71^^)i the fringe visibility must be 



V> 



_n_ 



(12) 



Time required to achieve equal precision: 
post-selected vs. non-post-selected 

In this section we calculate the number of experimental 
runs, which is equivalent to the time required to perform 
the experiment, which are necessary for a post-selected 
NOON state scheme to achieve the same level of precision 



as a non-post-selccted scheme when losses and imperfect 
detectors (characterized by rj and rj^) and other imper- 
fections (characterized by V and r]p) arc present. The 
apparent precision of a post-selected iVOOiV state scheme 
is given by 



A0P 



1 



ly'VN 



(13) 



where TV is the number of photons detected in coincidence 
events at the output of the interferometer. Instead, if we 
consider the number of photons N sent into the interfer- 
ometer as the relevant resource we obtain 



A(^„ 



1 



ufVN 



(14) 



By comparing Eqs. (fT3|l and ((T4)l it follows that we need 



(15) 



/ VpivVd) 



N 



more experimental trials in the post-selected scheme to 
obtain precision equal to the lossless case, which scales 
exponentially with the size of the A^OOTV-state. 

Similarly we can compare a iVOOiV-state scheme with 
a classical scheme of precision 



A(^SIL = 



1 



ViV^ 



(16) 



Comparing this with Eq. 
state scheme requires 



HI) we see that the NOON- 



(17) 



lyr] 



Vpivm)''v^N 



trials to achieve the classical precision. From this point of 
view the A^OOA^-state scheme beats the classical scheme 
only if rj / {rjp{rjrjtj)^ V''^ N) < 1, which is challenging to 
achieve in practice. 



Heralded Holland-Burnett state generation 

Holland-Burnett or twin-Fock states [l2| , which we de- 
note HB(7V = 2K), where K € {1, 2,3,.. .}, can be pre- 
pared by sending K photons to both inputs of a 50:50 
beam splitter. The state of the field at the two output 
modes of the beam splitter is given by 



|*k) = J2 A.mji\2m)a\2iK-m))t 



m=0 

where the coefficients are 



Am.K = (-1) 



K—rn 



/2m\ (2{K -m)\ /l 
V m j V K -m J \2 



2K 



(18) 



1/2 



This notation differs slightly from the standard notation, 
where K = N, and there are a total of iVtot = 2A^ pho- 
tons instead of N. We choose our notation to clarify 
the comparison between A^OOA^ states and the HB states. 
Note that the HB state has similar structure to the states 
examined for optimal phase estimation in lossy interfer- 
ometers in Ref. [J]. 

The heralded HB-state source is based upon generation 
of heralded Fock states by means of parametric down- 
conversion in Potassium Dihydrogen Phosphate (KDP) 
IJ] . A crucial feature of this nonlinear source is that the 
modes are entangled only in photon number and have 
no entanglement in any other degrees of freedom. This 
allows detection of N photons in one mode to herald, 
in principle, N photons in the other mode in a pure 
spatiotemporal quantum state. No spectral filtering is 
needed to increase the state purity of the heralded pho- 
tons. This scheme is therefore scalable, since there is no 
decrease in heralding efficiency arising from discarding 
photons by filtration. The use of high-efficiency number- 
resolving heralding detectors enables production of pre- 
cisely the target HB state abrogating the need for an- 
cillary assumptions about contamination from higher or 
lower photon numbers. 

A mode-locked Ti:Sapphire laser system operating at 
80 MHz provides 100 fs duration pulses centered at a 
wavelength of 830 nm. These pulses are frequency dou- 
bled in Beta Barium Borate (BBO) to a wavelength of 
415 nm. After filtering any residual infrared radiation, 
the UV-pulses are used to concurrently pump two KDP 
crystals oriented for type-II collinear parametric down- 
conversion (PDC) [14|. The orthogonally polarized PDC 
modes are split using polarizing beam splitters. The ver- 
tically polarized mode of each PDC source is detected 
by a fiber-coupled avalanche photodiode (APD) serving 
as a herald signal for the horizontal modes. After ro- 
tating one of the horizontal modes to vertical polariza- 
tion, the heralded modes are sent to a fiber-coupled po- 
larizing beam splitter (FPBS) where they are combined 
into a single spatial mode. Each heralded beam and the 
FPBS output are passed through polarization compen- 
sation wave plates to ensure appropriate polarization af- 
ter passing through the optical fibers. The FPBS out- 
put has additional polarization compensation. Success- 
ful preparation of the heralded state is signaled by a co- 
incidence detection event between these two heralding 
detectors. When only one photon is registered at each 
herald detector, this heralds the state |1,1)hv) where 
|7i,77t,)hv denotes n (m) horizontally (vertically) polar- 
ized photons. This is equivalent to the two-photon HB 
state (|2,0)+_ - |0,2)+_)/a/2, where -I-,- denotes the 
45° rotated basis. Note that single photons are typically 
heralded with an efficiency of 0.18 in this experimental 
configuration. 
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Interferometry 

We used a polarization-based interferometer to apply 
phase shifts between +45° and —45° plane polarization 
modes. This consists of a half-wave plate situated be- 
tween two quarter-wave plates. The quarter-wave plates 
were initially aligned to rotate + into — without the half- 
wave plate. The addition of the half-wave plate allows 
application of the phase, 0, being twice the angle of the 
half-wave plate. The + and — modes are then interfered 
on a polarizing beam splitter for analysis. 



photon number subspaces when measured in coincidence 
with the herald signal. The wave plate settings used, 
{0,,^}, were {0°,0°}, {0°, 11.25°}, {0°,22.5°}, {0°,45°}, 
{22.5°, 0°}, {22.5°, 22.5°}, {22.5°, 45°}, and {45°, 22.5°}, 
where 9 and are the angles of the QWP and HWP re- 
spectively. This forms an over-complete set of measure- 
ments on the whole space. The state was reconstructed 
by maximum likelihood estimation [2l| using the mea- 
surement outcomes. In order to check the consistency of 
the reconstructed density matrix, we calculate an overlap 
of 0.85 of the two-photon subspace with the ideal two- 
photon HB state. Fig. [St, in agreement with the Hong- 
Ou-Mandel interference visibility. 



State tomography of the heralded state 

We characterize the polarization and photon number 
degrees of freedom of the heralded state by performing 
state tomography [11, [lO]. A half- (HWP) is followed 
by a quarter- wave plate (QWP), the angles of which de- 
termine the measurement implemented. The HWP and 
QWP are followed by a polarizing beam splitter (PBS) 
with the transmitted mode coupled into a single mode 
fiber-coupled APD and the reflected mode coupled into a 
single-mode fiber beam splitter (50:50) with outputs con- 
nected to fiber-coupled APDs. This allows partial photon 
number resolution, giving information about the different 



Quantum Fisher information of the reconstructed 

state 



For the reconstructed density matrix in our experi- 
ment, the quantum Fisher information is given by [l9[ 



k,l 



{Pk - Plf 
Pk +Pl 



m\n\^i)? 



(19) 



where h is the photon number operator in the sensor 
arm of the interferometer, and {pk} and {|^fe)} are the 
eigenvalues and corresponding eigenvectors of the recon- 
structed density matrix. 



